Formation of Dark Matter Haloes in a Homogeneous Dark Energy Universe 
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Several independent cosmological tests have shown evidences that the energy density of the Uni- 
verse is dominated by a dark energy component, which cause the present accelerated expansion. The 
large scale structure formation can be used to probe dark energy models, and the mass function 
of dark matter haloes is one of the best statistical tools to perform this study. We present here 
a statistical analysis of mass functions of galaxies under a homogeneous dark energy model, pro- 
posed in the work of Percival (2005), using an observational flux- limited X-ray cluster survey, and 
CMB data from WMAP. We compare, in our analysis, the standard Press- Schechter (PS) approach 
(where a Gaussian distribution is used to describe the primordial density fluctuation field of the 
mass function), and the PL (Power Law) mass function (where we apply a nonextensive q-statistical 
distribution to the primordial density field). We conclude that the PS mass function cannot explain 
at the same time the X-ray and the CMB data (even at 99% confidence level), and the PS best 
fit dark energy equation of state parameter is cj = —0.58, which is distant from the cosmological 
constant case. The PL mass function provides better fits to the HIFLUGCS X-ray galaxy data 
and the CMB data; we also note that the uj parameter is very sensible to modifications in the PL 
free parameter, q, suggesting that the PL mass function could be a powerful tool to constrain dark 
energy models. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 

In the last few years, several experiments and obser- 
vations have shown strong evidences that the expansion 
of our Universe is accelerating, under the influence of a 
mysterious dark energy component [1, 2]. The dark en- 
ergy presents an equation of state parameter u = p/p and 
represents about 74% of the Universe. When uj = —1 we 
have the special case of the cosmological constant - or 
vacuum energy - as responsible for the Universe acceler- 
ation (it is the standard LCDM cosmological model). 

We also know that around ~ 26% of the Universe are 
composed by matter (baryonic + dark matter); analyzing 
the growth of these matter density fluctuation, we can 
describe the formation of the large structures, such as 
clusters and super-clusters of galaxies. 

The mass function of galaxies is a powerful tool to 
study the large scale structure formation in the Universe. 
The standard analytical approach was developed in 1974 
by Press and Schechter [3] (hereafter PS), and it is still 
used today, for its success and simplicity. The PS ap- 
proach uses a Gaussian distribution to describe the pri- 
mordial matter density fluctuation field. 

Inspired by Tsallis g'- nonextensive statistics Q and 
kinetic theory [E*, "g^, we have proposed a new mass 
function [7- 9], replacing the Gaussian by a non-Gaussian 
distribution - a Power Law (PL) distribution. The PL 
mass function has a physically motivated free q parame- 
ter (related with the long range gravitational correlations 
between particles), which provides malleability to fit ob- 
servational data. Also, if (7 = 1 we recover the original 



* Electronic address: 11uciomarassi@ect.ufrn.br] 



Gaussian distribution. The PL mass function is an ex- 
tension of the PS one, and presents the same analytical 
simplicity. 

We have compared the PS mass function with the 
PL one, in the special case of the LCDM cosmological 
model [9i]. Using the X-ray flux-limited sample of galaxy 
clusters from Reiprich and Boehringerploj (HIFLUGCS) 
in a statistical analysis, and applying independent 
cosmological tests (BAO and Shift Parameter) to bet- 
ter constrain the results, we concluded that the PS ap- 
proach presents incompatibilities with the independent 
CMB data. On the other hand, using the PL mass func- 
tion, we have better cosmological parameters, and we 
note an overlap with the CMB data for a large range of 
values of its q free parameter. These results, although 
encouraging, are limited to the LCDM model {uo = —1). 

The mass function could be used to probe the dark 
energy behavior in the Universe, and works using homo- 
geneous and non-homogeneous dark energy in the large 
scale structure formation are beginning to flourish in the 
literature. In this work we use a model of homogeneous 
dark energy from Percival By allowing other values 
for the UJ dark energy parameter, we can compare the ap- 
plicability of the PS and the PL mass functions to probe 
dark energy models. 



II. THE MASS FUNCTION EQUATIONS 

The distribution of bound objects with masses be- 
tween M and M+dM, using the Gaussian distribution to 
describe the primordial density fluctuation field 5 = 5p/ p 
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(The PS approach) reads Q: 
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where cr^^^j = ) is the mean squared fluctuation, and 
Sc is the critical S for collapse. 

Now, if instead of Gaussian we use the PL distri- 
bution for the initial fluctuations, we derive the follow 
expression j^-Q 
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where the factor Bq is a one-dimensional normalization 
constant. In the limit q ^ 1 the above PL expression 
reduces to the standard Gaussian approach. 



III. RESULTS AND FUTURE PERSPECTIVES 

We performed a statistical procedure using the X- 
ray HIFLUGCS data sample[10], and so we determined 
the statistical confidence contours (from la to 3cr errors) 
for the pair of parameters, and Vtrn- We applied the 
homogeneous dark energy model of Percival[ll], and we 
marginalized over all possible values of the dark energy 
equation of state parameter uj. 

In Figure [1] we show the contours in the VLj^-cf^ plane 
obtained by using the PS and the PL approaches. The 
left panel shows the PS results: the best-fit to the HI- 
FLUGCS X-ray galaxy cluster data are Vtra = 0.12, 
(78 = 0.98 and uo = —0.58. The best fit values of Qrn 
and (78 are, respectively, very low and very high when 
compared to the nowadays independent CMB results [2^. 
the solid vertical and horizontal lines show the minimum 
and maximum CMB limits (from WMAP 3 years data) 
for the Qrn and as parameters. The contours using the 
Gaussian distribution (the left panel) does not intercept 
the CMB values. In the other hand, by using the PL 
mass function with a certain q free parameter (the right 
panel of Figure [T]), we could overlap the WMAP indepen- 
dent data, for almost the same dark energy uj parameter 
(in this example, we use q = 1.15). In Figured] we see 
that, for almost the same particular model of dark en- 
ergy {uj ~ —0.5), we obtain better cosmological parame- 
ters using the PL mass function, while the standard PS 
mass function cannot explain the independent WMAP 
data even at 99% confidence level. 

In Fig. [21 we plot two statistical contours in the ^m-o's 
plane, using the HIFLUGCS data with the PL mass func- 
tion, under the dark energy model of Percival[ll]. Again 
we marginalize over all possible values of cu. By fixing 

= 1.10 in the PL mass function, we obtain as best-fit a 
dark energy uj = —0.48, Qrn = 0.15 and as = 0.85 (and 
the plot overlap the independent WMAP data with this 
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FIG. 1: Statistical contours fitting the mass function to the 
HIFLUGCS X-ray data. We plot the Qm—crg plane, marginal- 
izing over the dark energy cj parameter. In both panels, the 
solid vertical and horizontal lines show the minimum and 
maximum WMAP limits for Qm and as. Note that the con- 
tours using the Gaussian distribution (the left panel) does not 
intercept the WMAP values while the other one (in the right 
panel), using the PL distribution with q = 1.15, overlap the 
WMAP independent data, for almost the same dark energy 
cj parameter. 



configuration). When we fix g = 1.20 the plot continues 
to overlap the WMAP data, with almost the same best- 
fit for Qrn and (78, but it presents a huge modification in 
the best-fit dark energy u parameter {u = —1.08 here, 
close to the cosmological constant case). So, in the range 
1.10 < q < 1.20, the PL mass function provides better 
fits to the HIFLUGCS X-ray galaxy data and the inde- 
pendent data from CMB. Also, the sensibility over the cu 
parameter, as we change the free q parameter, indicates 
that the PL mass function could be a powerful tool to 
constrain dark energy models. 

It must to be noted, however, that we are using the 
standard Press & Schechter approach, in this present 
work. This approach smooths the initial density fluctua- 
tions, and uses the spherical model to calculate the epoch 
of critical overdensity for collapse of the density pertur- 
bations {6c - critical linear density contrast at collapse 
time). Such simple model eventually fail in detail, given 
structure formation complexities like, for example, the 
asymmetrical gravitational collapse; meanwhile, recent 
mass function improvements from numerical simulations 
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FIG. 2: Statistical contours on the Qm-crg plane, using the 
HIFLUGCS data with the PL mass function, under the dark 
energy model of Percival (2005). Marginalizing over the dark 
energy cj, we see that the PL mass function overlap the in- 
dependent CMB data (solid lines), for a wide range of its 
free parameter (1.10 < q < 1.20). Note that while q grows 
from 1.10 to 1.20, the parameter cj presents huge modifica- 
tions (from -0.48 to -1.08) - the cj parameter is highly sensible 
to changes in the PL mass function. 



have now quantified these problems - e.g. the Sheth & 
Tormen ellipsoidal collapse (1999 - hereafter ST)[l2|, and 
the universal function of Jenkins et al. (2001 - hereafter 
Jenkins) [13]. 

In this work, we use the Percival critical overdensity 
for collapse (Sc)^ with growth factor derived from the 
special case of flat dark energy models with constant uj 
parameter (equations 19 to 21 from reference ^3). This 
critical overdensity is only weakly dependent on cosmo- 
logical parameters, and the cumulative mass function, 
using the numerical fit of ST, presents very little differ- 
ence for different cosmologies (different uj dark energy 
parameters) at low redshifts {z close to zero); but as we 
go further back in time (towards higher redshifts), the 
difference becomes appreciable (figure 6 from [11]). We 
see that the evolution of the mass function is strongly de- 
pendent on cj, due to the evolution of the linear growth 
factor (which affects Sc directly). 

It is worth to discuss, at this point, the results from the 
work of Matthew et al. (2009 - hereafter Matthew) [14]. 
Matthew studied the role of the linear density contrast at 
collapse time, (5c, in Early Dark Energy (EDE) models - 



where dark energy has a non-negligible importance since 
the beginning of the structure formation process. Us- 
ing a pure numerical approach to compute (5c, Matthew 
found that, at redshift z = 0, the EDE mass functions 
are not greatly altered compared to the LCDM cosmo- 
logical model, but the difference increases with redshift 
- the same conclusion as Francis et al. (2008, hereafter 
Francis) [15]. Matthew also demonstrates that the EDE 
model presents a basic agreement between the ST and 
the Jenkins mass functions, even if we use a cosmology 
dependent (5c, derived from the spherical collapse model. 
However, when we compare the Jenkins mass function 
(which is 'blind' to the growth history of the universe) 
to the ST mass function with a fixed 5c {Sc = 1.689), 
or deriving ST Sc from the spherical model (where the 
history of the universe - the growth history - counts), we 
can see small differences between these mass functions in 
the high mass end (M > 10^^ M^/h); despite these dif- 
ferences, the most up-to-date numerical simulations have 
not sufficient accuracy at this high mass range to dis- 
criminate between these two approaches. 

So, we observe a general agreement between the con- 
clusions of Percival, Matthew and Francis, concerning 
dark energy models in the structure formation. All dark 
energy mass functions, ST or Jenkins, have almost the 
same behavior as the LCDM model at z = 0, and at high 
redshifts we can perceive differences between the cosmo- 
logical models, due to the evolution of the linear growth 
function. In the work of Percival, the focus is the dif- 
ference based on the dark energy oo parameter, and in 
the work of Matthew and Francis, the focus is the differ- 
ence between the ST (with fixed or evolving Sc) and the 
Jenkins mass function. 

In this present work we use, as said before, the stan- 
dard PS approach, with a Gaussian and a non-Gaussian 
(PL) initial distribution function. We use in this work 
a low-redshift X-ray galaxy cluster survey, but even at 
small redshifts we show a strong mass function depen- 
dence on the dark energy uj parameter, using the PL 
approach, and the statistical procedure show signif- 
icant contour differences between the dark energy mass 
function and the LCDM cosmological model. This show 
clearly the power of the PL mass function q free param- 
eter to fit observational data; and if this malleability is 
observed to galaxy clusters at 2: ~ 0, we should expect 
even more exciting results using high redshifts catalogs, 
to probe the characteristics of the dark energy, in a close 
future. Also, we intend to study these dark energy mod- 
els using ST and Jenkins mass functions as well, studying 
the Sc influence in the process, to compare our results 
with the previous conclusions of Percival, Matthew and 
Francis. 



IV. CONCLUSIONS 

In this work we perform a statistical analysis, using the 
X-ray galaxy cluster data named HIFLUGCS [10], and 
the independent Cosmic Microwave Background radia- 
tion (CMB) data from WMAP[2], to study the influence 
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of galaxy mass functions under a homogeneous dark en- 
ergy background, from the work of Percival[ll] . Perform- 
ing a analysis and marginalizing over the uo = p/ p dark 
energy equation of state parameter, we compare the stan- 
dard Gaussian Press-Schechter mass function [3] to the 
PL mass function (based on the non-extensive thermody- 
namics and kinetic theory), which presents a physically 
motivated free parameter, g0-[2|- 

We conclude that the PS mass function cannot explain 
at the same time the X-ray and the CMB data, even at 
3(7 confidence level, which imposes a problem to the stan- 
dard Gaussian scenario. Also, the PS approach presents 
a best fit dark energy equation of state parameter of 
uj = —0.58, which is clearly distant from independent 
estimates of a; — 1 (a cosmological constant). 

Observing the contours of Figure[2l in the range 1.10 < 
q < 1.20, the PL mass function provides better fits to 



the HIFLUGCS and the CMB independent data, with 
almost the same best-fit for and ag, but with a huge 
modification in the best-fit dark energy cu parameter, as 
we change q. This strongly suggests that the PL mass 
function could be a powerful tool to constrain dark en- 
ergy models. Future analysis, based on the next genera- 
tion of cluster surveys in medium and high redshifts, will 
certainly provide better constraints to the dark energy 
model that accelerates de Universe. 



Acknowledgements 

The author are grateful to ECT-UFRN and FAPESP, 
for the financial support provided, and to Lima J.A.S. 
for helpful discussions. 



[1] A. G. Riess et a/., Astrophys. J. 659, 98 (2007); R. G. 
Carlberg et al, ApJ. 462, 32 (1996); A. Friaca, J. S. 
Alcaniz and J. A. S. Lima, Mon. Not. R. Astron. Soc. 
362, 1295 (2005), ArXiv:0504031 [astro-ph]; S. W. Allen, 
R. W. Schmidt and A. C. Fabian, MNRAS 334, Lll 
(2002); J. A. S. Lima, J. V. Cunha and J. S. Alcaniz, 
Phys. Rev. D 68, 023510 (2003), ArXiv:0303388 [astro- 
ph]; J. V. Cunha, L. Marassi and R. C. Santos, IJMPD 
16, 403 (2007); J. V. Cunha, L. Marassi and J. A. S. 
Lima, MNRAS 379, L1-L5 (2007), ArXiv:0611934 [astro- 
ph]. 

[2] D. N. Spergel et al., Astrop. J. Suppl. 148, 175 (2003); D. 

N. Spergel et al., Astrophys. J. Supl. 170, 377 (2007); J. 

Dunkley et al., ArXiv:08030586 [astro-ph]; E. Komatsu 

et a/., ArXiv:08030547 [astro-ph]. 
[3] W. H. Press, P. Schechter, ApJ 187, 425 (1974). 
[4] C. Tsallis, J. Stat. Phys. 52, 479 (1988). 
[5] R. Silva, J. A. S. Lima and J. Santos, Phys. Rev. E61, 

3260 (2000); F. Valentini, Phys. Plasmas 12, 072106 



(2005). 

[6] J. R. Bezerra, R. Silva and J. A. S. Lima, Physica A 322, 
256 (2003). 

[7] J. A. S. Lima and L. Marassi, IJMP D13, 1345 (2004). 
[8] L. Marassi and J. A. S. Lima, IJMP D16, 445 (2007). 
[9] L. Marassi, J. V. Cunha and J. A. S. Lima (2008), 

ArXiv:08070444 [astro-ph]. 
[10] T. H. Reiprich and H. Boehringer, Astrophys. J 567, 716 

(2002). 

[11] W. Percival, A&A, 309, 823 (2005). 

[12] R. K. Sheth and G. Tormen, MNRAS, 308, 119 (1999). 

[13] A. Jenkins, C. S. Frenk, S. D. M. White, J. M. Colberg, S. 

Cole, A. E. Evrard, H. M. P. Couchman, and N. Yoshida, 

MNRAS, 321, 372 (2001). 
[14] J. F. Matthew, F. L. Geraint and V. L. Eric, MNRAS, 

393, L31 (2009). 
[15] M. J. Francis, G. F. Lewis and E. V. Linder 

(2008),ArXiv:08082840 [astro-ph].. 



